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Abstract 

The general formalism of the multipolar expansion of electrostatic interactions is applied to the 
calculation of the potential energy between a excited atom (without fine structure) and a ground 
state diatomic molecule at large mutual separations. Both partners exhibit a permanent quadrupole 
moment, so that their mutual long-range interaction is dominated by a quadrupole-quadrupole 
term, which is attractive enough to bind trimers. Numerical results are given for an excited 
Cs(6P) atom and a ground state Cs2 molecule. The prospects for achieving photoassociation of a 
cold atom/dimer pair are thus discussed and found promising. The formalism can be generalized 
to the long-range interaction between molecules to investigate the formation of cold tetramers. 
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I. INTRODUCTION 



Since it was proposed by Thorsheim et al. l| in 1987, and first observed for Sodium 
and for Rubidium 3| atoms in 1993, tlie photoassociation (PA) of pairs of ultracold atoms 
has had a tremendous impact on research in atomic, molecular, and optical physics at low 
teinperatures. There are several recent review articles devoted to the various aspects of PA 
#iZ|, therefore we briefly recall below some of the main features of the PA process, which 
gave rise to a new high-resolution spectroscopic technique, i.e. the PA spectroscopy. Due to 
their extremely low relative kinetic energy, atoms from an ultracold gas can be associated 
via a quasi-resonant free-bound dipolar transition to form an electronically excited molecule, 
which is often created in a highly excited rovibrational level. As the PA process is mainly 
controlled by the long-range electrostatic interactions between cold atoms, it has been used as 
a high-resolution spectroscopy technique for highly rovibrational levels. The highly excited 
rovibrational levels observed using PA correspond to vibrational motion of a molecule with 
much larger extension than the usual chemical bond 8l-ll2l|. Such molecules with a very large 
amplitude of vibration had been predicted fifteen years before the mentioned experiments 
HQ . The spectroscopy of the highly excited rovibrational levels of photoassociated dimers 
permitted, in particular, to determine the most accurate values of the radiative lifetime of 



the first excited state of alkali- metal atoms (see for instance 15|]). Another example of PA 



application is the formation of stable ultracold molecules, reported initially for Cs2 16 [, and 



17|-|24|. 



later for many other homonuclear and heteronuclear alkali-metal diatomic molecules 

With the improvement of the experimental techniques at ultra-cold temperatures, the 
study of the quantum dynamics of few-body systems in the ultracold regime has become pos- 



sib 



e, as illustrated by the recent observations of cold collisions between atoms and molecules 



25N29| . Such phenomena attract at present a lot of interest as they represent the first man- 



ifestation of a novel ultracold chemistry, which is controlled by the quantum nature of the 
colliding partners 30|, |31|. In particular, at certain conditions, the ultracold few-body dy- 
namics exhibits universal (i.e. species-independent) properties for long-range bound states 
and resonances (see, for example, Ref. 
as the Efimov physics 



34 



and references therein), nowadays referred to 
351]. The Efimov states have recently been observed experimentally 



All these developments concern atoms and molecules in their electronic ground state. 



The purpose of the present study, as the first of a series of papers, is to investigate the next 
step towards ultracold chemistry: the association of ultracold atoms and molecules with a 
laser field to create weakly bound trimers or tetramers in an excited electronic state, which 
has not been previously discussed in the literature to our knowledge. Just like for pairs of 
atoms, the PA probability is determined by the long-range interactions between the colliding 
partners. Here, we consider the long-range interaction between a molecule in a given 
rovibrational level {vd,j) with an atom in a P electronic level without fine structure. This 
situation will be illustrated with the interaction between a ground state Cs2 molecule, and 
an excited Cs(6P) atom. The leading term of this interaction at large interparticle distances 
i? is a quadrupole-quadrupole term varying as R~^, but can be easily generalized to other 
species. The present work can also be viewed as a step beyond several related studies. The 
quadrupole-quadrupole interaction between two exited atoms has been calculated for 



alkali-metal atom pairs 39| and for the LiB molecule 40|. In Refs. 4l|, |42|, the van der 



Waals interaction (varying as R ^) between alkali-metal dimers in the {vd = 0,j = 0) level 



their lowest 
in Refs. 43 



triplet state and a ground state alkali-metal atom has been determined, while 
44| the interaction between a ^11 molecule and a atom at fixed geometries is 



obtained as a sum of a dipole-quadrupole term (in R ^) and a quadrupole-quadrupole term 
(in R-^). 

In Section ini we briefly review the main ingredients of the perturbative approach based on 
the multipolar expansion of the long-range interaction between the two fragments. Section 
mil is devoted to the calculation of C5 coefficients of the long-range behavior of molecular 
potentials. We consider the general case of an arbitrary rotational state j of the dimer as 
well as we give an analytical solution for the particular case of j = 1. Atomic units (a.u.) 
for distances (1 a.u. = 0.0529177 nm) and for energies (1 a.u. = 219474.63137 cm~^) will 
be used throughout the paper, except otherwise stated. 



II. INTERACTION POTENTIAL AND PERTURBATION THEORY 



We start the description of the present theory from the general for instance in 



Refs. 



43 



45| . We consider two charge distributions, A and B, far from each other such 



that they do not overlap with each other. A criterion for such a condition is given by the 
so-called Le Roy radius Q, ^ defined as Rlr = 2{^y {r\) + ^/{rB)), where {r\) and (r^) 



are the averaged squared distance of the outermost electron from the origin of each charge 
distribution A and B, respectively. The electrostatic potential energy of interaction between 
A and B can be written as an expansion over products of multipole moments of A and B 
located at a distance R from each other 

+00 L< ^ 

La,Lb=0 M=~L< 
X fLALBMQLA^A)Q~LB(^B), (1) 

where L< = min(Lyi, Lb)- The operator Qi^^{rx) is associated to the 2^^-pole of the charge 
distribution X {X = A or B), expressed in the body-fixed coordinate system with the origin 
at the center of mass of X 



Qf. irx) = J E g.rf-y£ 0,) , (2) 

where qi is the value of each charge i composing the distribution X. The two coordinates 
systems (centered at A and B) are assumed to have parallel axes with the Z axis that goes 
from the center of mass of A towards B (see Fig. [1]). This choice of Z implies in Eqs. ([1]) 
and ([2]) that Ma = —Mb = M, where Ma and Mb are the projections of La and Lb, so 
that the factor /laLbM is equal to 

(-1)^" (La + LbV. 

^"""^'^ y/{LA + M)\{LA-M)\ 

X ^ . (3) 

^iLB + M)\{LB-M)\ 

The energy of interaction between the two charge distributions is calculated using per- 
turbation theory. To the lowest (zeroth-) order of perturbation theory, the two systems are 
independent and the total energy is the sum of the individual energies 

-^0 = -^AO + ^BO ) (4) 
and the total wave function is the product of individual wave functions 

\<) = Ko)Ko)- (5) 

In Eqs. @, ([5]) and below, the superscript labels the perturbation order, and the subscript 
labels the unperturbed states. 



In the present study the system A is the alkah-metal dimer and B is the alkah-metal 
atom. We consider the dimer in a vibrational level Vd of its fundamental electronic state 
|X^E+,f(i), and in an arbitrary rotational state |j, mj). In order to investigate a realistic 
approach for atom-molecule photoassociation, we consider the atom B with a single outer 
electron being excited to the p state \n,i = 1, A). However, we ignore in the following the 
fine structure of the excited atom for clarity, as discussed later in the text. The projections 
rrij and A are defined with respect to the Z axis. The energy origin corresponds to an infinite 
separation between the atom and the dimer. Thus, the unperturbed energy reduces to 

E° = B^j (j + 1) , (6) 

where By^ is the rotational constant of the dimer in its vibrational level Vd- The atomic 
state of B is expressed in the LS coupling case, because the operators in the interaction 
potential of Eq. act only on the coordinate part of wave functions. The first order 
correction Eq to the energy is due to the permanent multipoles of A and B. In our case, 
both distributions exhibit a permanent quadrupole moment in their body-fixed frame, so 
that the most important contribution comes from the quadrupole-quadrupole interaction 
with an asymptotic coefficient C5 



III. CALCULATION OF THE C5 COEFFICIENT 

The C5 coefficient is calculated for arbitrary values of rrij and A, using the degenerate 
perturbation theory. We define two body-fixed coordinate systems (CS) (Fig. [1]). The first 
CS (we call it the dimer CS, or D-CS) with axes X^, Y^, and has as the origin the center 
of mass of the dimer. The axis is the dimer axis and the axis is orthogonal to the 
plane of the trimer. The second CS (trimer CS, or T-CS) with axes XYZ is such that the 
X, Z axes are also (as Xa and Za) in the plane of the trimer, while Z is oriented from the 
center of the dimer towards the atom B; the Ya and Y axes are identical. The T-CS is 
deduced from the D-CS by a rotation with an angle S around the Y axis. 

The perturbation Hamiltonian V^^(-R) for the quadrupole-quadrupole interaction is given 

5 





FIG. 1. The two coordinate systems, X^YaZa (D-CS) and XYZ (T-CS) defined for tlie dimer 
and for the trimer, respectively. The Y and axes coincide and point into the plane of the 
figure. The subsystem A in this figure is the Cs2 molecule, the subsystem B is the Cs atom. The 
T-CS is related to the laboratory coordinate system (xyz) by the usual Euler angles (a, /3, 7), not 
represented here. 



by setting La = Lb = 2 in Eq. ([T]): 



ab{ ) ^5 2^ (2 + M)!(2-M)! 



(8) 



M=-2 



The Hamiltonian has the form of a sum of tensor products, composed of operators acting in 
subspaces of the unperturbed eigenstates of A and B. Usually, the 2'^-pole tensor compo- 
nents qff' of a charge distribution are defined in its proper CS, i.e. XaYaZa for the dimer 
A. Therefore, the QfJ tensor components in the T-CS are written as 

Q¥= E dUAS)qf, (9) 



M': 



where dpjpp{S) are the reduced Wigner matrix elements. In the case of the alkali- metal 
dimer in the state, the only non-zero component of the quadrupole moment is ^2 



Eq. (Q reduces to 

Q2=dumi (10) 

The component is just a scalar parameter, which will below be referred to as 

In the T-CS, the wave function of the rotational state of the dimer \jmj) is written as 



A/(2j + 1) /2dl^^Q{6), depending only on the internal angle 6. The normalization constant 
is such that the integral over angle 6 is unity. We obtain the following expression for matrix 
elements of the operator (^2^ 



2j 



-4 



(11) 



where the Clebsch-Gordan coefficients Cf^l'^t appear after intergrating the product of three 



^MM' functions 



48 



The zeroth-order energy Eq depends on j and is degenerate for all 
values of mj. Thus, the perturbation Hamiltonian (Eq.®) has to be evaluated with the 
degenerate perturbation theory. Indeed, Eq. ffTTj) shows that the quadrupole moment has 



matrix elements for different values of mj because of the m'- = rrij + M selection rule. The 
degeneracy between different mj values will be removed, leading to different values of C5. 
This is the key point of the present treatment, as anisotropic values of C5 are determined 
as functions of quantum numbers of the partners, and not restricted to a given geometry. 

Assuming the alkali-metal atom being in the state labeled \niX) (the spin is neglected 
here), we calculate the matrix elements of the atomic quadrupole moment operator for a 
given i between two different Zeeman sublevels A and A' following the same treatment as 
above. We obtain 



(niX' 



Q 



M 



nlX^ 



Atx 




2tt 



(rle) / # 



(12) 



where the negative sign comes from the electron charge. The mean squared position (r,^^) 
of the valence electron is independent on A. Using the properties of spherical harmonics, we 
rewrite Eq. (fT2|) as 



(^n£X' 



Q 



niX 



^2Qm^2MlX \' nil ■ 



(13) 



The situation is analogous to the molecular case: if M 7^ 0, the operator couples A to 
A' = A + M, and the perturbation Hamiltonian lifts the degeneracy with respect to A also. 

Summarizing the above results, the perturbation operator of Eq. (|8]) couples the (2j + 1) 
rotational states of the molecule with a given value of j, and the (2£ + 1) Zeeman states of 



the atom with a given value of The C5 coefficients are then given by (2j + 1) x [2i + 1) 
eigenvalues of the operator V^%. Using Eqs. f lTTj) and f|T3|) . the matrix elements of V^^ are 
written 

{jm'/\'\VZ\jm,i\) = -24C^o°,A^^ 

\ ^ '^2Mjmj'^2-Me\ 

X 



(2 + M)!(2-M)!- 

From the integration over Euler angles and the properties of Clebsch-Gordan coefficients, 
the following selection rules for are derived: (1) The projection rhj oi the total orbital 
momentum J = j + £ on the laboratory z axis is conserved. (2) The projection m j = rrij + X 
of the total orbital momentum J on the Z axis of T-CS is conserved. This rule can also be 
deduced by the combination of Eq. ( fTTi) and Eq. ( |T3l) . 

Equation (HM demonstrates the equivalence between the atomic orbital momentum i and 
the dimer rotation j in the formalism, which describes the long-range interaction between two 
charge distributions with defined angular momenta irrespective to their internal structure. 
If one of the two angular momenta is zero, the corresponding quadrupole moment vanishes, 
and the C5 coefficient as well. Therefore, the interaction will be the usual Cq/R^ van der 
Waals potential. If neither of the two angular momenta, j and i is zero, the long-range 
interaction varies as C^/R^ and, therefore, the potential has a larger density of vibrational 
states close to the dissociation limit than the lowest electronic state of the system when 
the atom in its ground S state. Such a situation is favorable for the photoassociation of 
atom-molecule pairs into excited trimers, just like for the photoassociation of identical atom 



pairs (see for example, experimental work of Ref. 49|). 



IV. RESULTS AND DISCUSSION 



To illustrate the previous formalism, we ffist consider analytically the simplest case j = 
i = 1. The perturbation Hamiltonian reduces to a 9 x 9 matrix with elements calculated 
from tensor products of the atomic and dimer states. For simplicity, we omit the j and £ 
labels in the following, and the quantum states of the atom-molecule pair are denoted by 
projections {|mj. A)} only. All such states form the basis of the representation. If we sort 
the states by values of the conserved projection of the angular momentum rrij = rrij + A, we 
obtain the matrix of in a block- diagonal form. 

8 



The two blocks defined by Irrij, A) = |— 1, — 1) and |+1, +1) (mj = ±2) reduce to a single 
element with a negative value of the corresponding coefficient 

= ^ . (15) 

It produces an attractive interaction. Two other 2x2 blocks (with rrij = ±1) are defined by 
the two subspaces {| — 1,0) ; |0, —1)} and {|0, 1) ; |1,0)}. The corresponding C5 coefficients 

24(72 (^^ ) 

are — " (positive value) and zero. Finally, the last 3x3 block comes from the subspace 
{| — 1, 1) ; |0,0) ; |1, —1)} {jnj = 0). Two of the corresponding C5 coefficients are zero, and 
the third one is 

C5 , (16) 

with the eigenvector ^ (|-1, 1) + 2 |0, 0) + |1, -1)). The coefficient in Eq. ^ is negative 
with the largest magnitude out of all C5 coefficients obtained in the case oi j = i = 1. It 
corresponds to the most attractive configuration between the atom and the dimer, and is 
expected to be the most favorable for the photoassociation. 

The results of the calculation for the case j = £ = 1 are summarized in Table HI The 
second column of the table gives the eigenvectors of the Hamiltonian of Eq. ([8]) in 
the J = 1 subspace. The third column gives the so-called reduced values of C5 in units 
of 52 ('"nf)' which stress the general character of our treatment: it can be applied to all 
alkali-metal trimers, but it can also be compared with the existing results on the long-range 
interaction between two excited atoms |39|. The eigenvectors obtained here are the same as 
in Ref. 39|, but the signs of the C5 coefficients are opposite to the coefficients obtained in 
Ref. {391. The reason is clear from Eqs. ( fTT]) and ( fT2|) : The signs of the matrix elements of 
quadrupole moments for the dimer Eq. (11 II) and the atom Eq. (fT2il are opposite. When they 
are combined together in Eq. (JHl) they give an additional negative sign to the perturbation. 
The two matrix elements of the atomic quadrupoles give the positive sign to the perturbation 
matrix elements. 

The fourth column in Table |T] displays estimates for the C5 coefficients for Cs2+Cs. To 
the best of our knowledge, there are no available values for the quadrupole moment of Cs2 in 
the literature. Therefore, we calculated it for the electronic ground state with the Gaussian 



quantum chemistry package SOj using the MP2 method with the Def2-TZVPP basis 51 |. 



To check the accuracy of such an estimation, we first calculated the quadrupole moment of 



K2 and compared it to available accurate ab initio calculations 

9 



52|. We obtained 12.258 



mj 



-2 
-1 
-1 




+1 
+1 

+2 



1-1,-1) 
^(|-l,0) + jO,-l)) 
^(1-1,0) -|0,-1)) 
^(|-1,1)+2|0,0) + |1,-1)) 
^(|-1,1)-|0,0) + |1,-1)) 

;fe(|-l,l)-|l,-l)) 
^(|1,0)-|0,1)) 
75(|1,0) + |0,1)) 

ll,l) 



{<i {rid) 



25 

24 
25 



36 
25 



24 
25 



_6_ 

25 



Cs (a.u.) 



-279 
1116 


-1674 



1116 


-279 



TABLE I. Values of the C5 coefficient and their corresponding eigenvectors characterized by their 
mj vahie, for Cs2{X^T,^ ,Vd = 0,j = l)-|-Cs(6P). The values of C5 are given in units of q2 {f^f) 
in the second column, and in atomic units for Cs2+Cs in the third column. For cesium, the data 
are: {r^p) = 62.65 a.u. and (/g = 18.56 a.u. (see text). Due to the uncertainty over the results 
are given with a precision of 1 a.u. 



a.u., which differs by a factor 1.28 from the value 15.689 a.u. of Ref. 52|. For cesium, the 



Def2-TZVPP basis jsil contains also effective core potentials (ECPs) standing for the 54 
inner electrons of the core. We obtained for Cs2 the value of 14.51 a.u. that we multiplied 
by the same factor to estimate the Cs2 quadrupole moment to = 18.58 a.u. The mean 
squared radius of the 6P orbital of cesium, which is 62.65 a.u., is calculated using a Dirac- 
Fock method Issl- It is worth to mention that the values of (75_ shown in the table are of the 



same order of magnitude as the values for Cs(6P)-|-Cs(6P) 39 1. 

For £ = 1 and arbitrary j, the perturbation Hamiltonian is a 3(2j -|- 1) x 3(2j -|- 1) matrix, 
that can be diagonalized numerically. The eigenvalues obtained numerically for j = 2 to 4 
are given in Table HIl The C5 coefficients are of the same order of magnitude as for j = 1, 
but in average they become smaller in magnitude as j increases, due to smaller Clebsch- 
Gordan coefficients. The C5 coefficients are sorted by values of which, in analogy to 
diatomic molecules, are labeled S, 11, A, $, F, and H for \mj\ = to 5. For S states, the 
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symmetry 


j 


Cs (a.u.) 


symmetry 


J 


Cs (a.u.) 


S+ 


2 


-913 


A 


2 


-140 




2 


116 




2 


1136 




3 


-796 




3 


-835 




3 


145 




3 


-87 




4 


-755 




3 


736 




4 


157 




4 


-721 




2 


399 




4 


-11 




3 


465 




4 


623 




4 


489 




2 


-399 


n 


2 


-964 




3 


-245 




2 


-19 




3 


1175 




2 


584 




4 


-783 




3 


-783 




4 


-161 




3 


64 




4 


835 




3 


532 


r 


3 


-465 




4 


-739 




4 


-320 




4 


108 




4 


1208 




4 


522 


H 


4 


-507 



TABLE II. The C5 coefficients of the Cs2{X^T,^ ,Vd = 0, j)-|-Cs(6P) long-range interaction calcu- 
lated numerically for j = 2 to 4. C5 are sorted by projections mj = rrij -|- A of the total orbital 
angular momentum on the Z axis, and by the sign -|-/- of the wave function with respect to a 
reflection through the plane containing the Z axis. In analogy to a diatomic molecule, the eigen- 
states are labeled with S"*"/", 11, A, <I>, T, and H for mj = 0, 1, 2, 3, 4, 5, respectively. The values for 
(r|p) and ^"^^ the same as in Table [H 

reflection symmetry through the Z axis is also considered, giving the usual +/- superscripts. 
For states other than S the sign -|-/- is not specified because such states are degenerate (in 
the present approximation) with respect to the reflection. 

We use the same symmetry notations in Figs. [21 [3] and HI where we display the long- 
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range potential energy curves C^/R^ for the Cs2+Cs(6P) system calculated for the first 
five rotational levels j of Cs2 as a function of the atom-dimer distance R. The energies of 
dissociation are given by Cs2 rotational energies, Bqj (j + 1), j = 0, . . . , 4. The rotational 
constant for ground vibrational level of Cs2 is Bq = 1.17314 x 10~^ cm~^ 541 ] . 

The potential energy curves are shown up to /2 = 500 a.u. Beyond this limit, the 
distance between Cs2 and Cs becomes comparable to the wavelengths of relevant atomic 
and molecular transitions, which are in the optical frequency domain. Therefore, in that 
region, electrodynamics effects, for example, due to retardation, should be taken into account 



55| 



As already mentioned, the lower limit of the region where the present approximation is 
applicable can be estimated by the Le Roy radius Rlr = 2 ^a/ (rQ(Cs2)) + a/ {i^lpiCs))^ , 
where {r^) and {r^p) are related the extension of the dimer and atomic electronic clouds, 
respectively. For the atom, one has (rip) = 62.65 a.u., which is given in Table [B As for the 
dimer, (rg) is calculated from the elements of the quadrupole tensor. All its non-diagonal 
elements Qa/3 are zero for states in D-CS. Its diagonal elements Qaa (a = ^AyYA, or 
Za) in D-CS are estimated using the Gaussian package, and is formally written as a sum 
over all charges 

Qaa = ^qi{a'f)- (17) 

i 

As the two nuclei of Cs2 are along the axis and it is a molecular state, QxaXa 
is equal to QyaYa both are functions of the coordinates of two valence electrons only 
(^ = 1,2) 

2 

Qx,x, = -eJ](Xi> . (18) 

1=1 

Now considering for simplicity that the cores are fixed at the position Z^ = i^^e/S (the rigid 
rotor approximation, valid for Vd = 0), for QzaZa obtain 

QzAZA = e"-^-ef2{Zl). (19) 

1=1 

Now, setting 

(-o^) = E E (20) 

i=l a={XA,YA,ZA) 

we reach the final expression for (r^) 

2 

(rl) = y - QzaZa - 2QxaXa , (21) 
12 
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FIG. 2. Long-range potential energy curves C^/B? as a function of the atom-dimer distance R 
(notice logarithmic scale along R), for the and symmetries, and for the five lowest rotational 
levels of Cs2(X^S+('yrf = 0))- The curves are drawn for distances larger than the Le Roy radius 
Rlr = 45 a.u. 

where e = 1 in atomic units. The Q^a niatrix elements are calculated with the same 
method as ^2 (using the ratio 1.27 to the K2 value) [56|, which yields QxaXa = "~69 a.u., 
QzaZa = —41 a.u. and = 8.7 a.u. Therefore, we obtain {r^) = 216 a.u. and the Le Roy 
radius Rlr = 45 a.u. 

As we can see from Figs. [21 13] and HI the Le Roy radius is smaller than the distance at 
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which the curves start to cross each other. This is the second main result of the paper. Unhke 
the case of two atoms, the rotational structure of the dimer is small enough to compete with 
the quadrupole-quadrupole interaction. The lower limit of for the applicability of the 
present perturbation approach is thus fixed by the crossing of the potential energy curves. In 
order to estimate Rm, we note that the first crossing occurs between the curves dissociating 
towards to the j = limit and from the most attractive curve corresponding to j = 1. 
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FIG. 4. Same as Fig. [2] for $, P and H symmetries. The upper panel displays $ symmetry, and 
the lower one displays both F symmetry (solid lines) and H symmetry (dashed lines). 

Putting 2Bo = C^/Rf^ yields a general estimate for Rm (see Eq. (ITG!) ) 



Rm ~ 



2Bo 
0.936 X 



25^0 

.0 W 1/5 



5, 







(22) 



For cesium, Eq.f l22p yields R ~ 102 a.u. In the rigid rotor approximation with Bq = 1/ (2/irg) 
and where n is the reduced mass of the dimer, Eq. fl22|) shows that the value of Rm is smaller 
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for lighter atoms. For example, replacing; in our treatment Cs by ^Li with the atomic 



parameters (r|p) = 32.5 a.u. |57|, = 5.05 a.u. and ^2 = ^^■'^ ^-^ [521, we obtain 



Rm = 43 a.u. This value of Rm is larger than the Le Roy radius for lithium, for which we 
obtained 26 a.u. using Eq. fl^T]) . 

For distances such that Rifj < R < Rm., the long-range potential ([T]) is still valid, but 
not the perturbative approach. The non-adiabatic interaction at the curve crossings (for a 
given symmetry) is expected to be strong. In particular, the interaction between permanent 
quadrupoles would couple the dimer rotational level j with j' = j±2, j±4 near the crossings. 
Higher-order contributions in should also be considered. 

The number of partial waves involved in the atom-dimer collisions depends on the 
temperature in the actual experiment. In order to give an upper bound for A^, we consider a 
potential curve with the most attractive C5, given by Eq. f|T6|) and with the added centrifugal 
term. It is straightforward to show that the height of the potential barrier E^""^ for a given 
N is 

^f5_Y' / iv(iv+i) y/^ 



^0.15Sx(^)"(.0K.»--. (23) 

Converted to the temperature, E^""^ is approximately 1 fiK for cesium. If we take typical 
temperatures 10-100 fiK for which photoassociation experiments are achieved, only a few 
partial waves (6~7 for the present case) will play a significant role in the collision. This 
contrasts with the PA of identical atom pairs, interacting with a long-range R~^ potential 
which allow much more partial waves than here. 



V. CONCLUSIONS AND PERSPECTIVES 



In this article, we used the multipolar expansion to calculate the long-range interaction 
energy of a diatomic molecule in its electronic ground state and in an arbitrary rovibra- 
tional level, and an excited atom. We applied our treatment to the case of a ground state 
Cs2 molecule and an excited Cs(6P) atom, as a prospect for cold atom/molecule photoas- 
sociation. The dimer and the atom interact through their permanent quadrupole moment. 
In contrast with previous works, the anisotropic interaction is computed for arbitrary ge- 
ometries of the atom/molecule pair, and depends on their internal quantum numbers. We 
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showed that the interaction hfts the degeneracy over their respective magnetic sublevels. 
Using the degenerate perturbation theory, we calculated the C5 coefficients characterizing 
the quadrupole-quadrupole interaction, for the five lowest rotational levels of the ground 
state dimer. The photoassociation of a ground state ^^S^ alkali-metal dimer molecule with 
a ground state nSi/2 alkali-metal atom is found possible by exciting the dimer-atom system 
with a laser frequency red-detuned from the nS — )■ nP atomic transitions. 

We demonstrated that the small- i? limit of applicability of our treatment is not due to the 
overlap of the electronic clouds of the partners as in the atom-atom case, but to the competi- 
tion between the rotational energy of the dimer with the long-range quadrupole-quadrupole 
interaction. This induces crossings between potential energy curves corresponding to differ- 
ent rotational levels. In the region 50 < i? < 100 a.u., the multipolar expansion is still valid, 
but not the perturbation approach. The inclusion of non-adiabatic couplings is required in 
this region for an appropriate description of the long-range behavior, as well as higher-order 
effects in 1/ R. This will be discussed ind detail a forthcoming paper. 

It is important to stress that the above treatment has been developed in the framework of 
the LS coupling case, in order to keep our description simple. The next step is to account for 
the fine structure of the excited atom. The main difference with the formulas reported here 
will be the change of the atomic state from P = \n,i, A) to Pj = \n,i,j, Xj) (with j = 1/2 
or 3/2), where the state Pj is written as the appropriate superposition of atomic states with 
different A and spin projections through a unitary transformation. We note however that 
for most of the alkali-metal atoms (from Na to Cs) the fine structure is much larger than the 
magnitude of the long-range atom-dimer interaction. Therefore, the related C5 coefficients 
will result from linear combinations of the coefficients of Table [Tll and will not modify the 
main statement of our study concerning the range of validity of our approach. In contrast, 
the case of Lithium atom will be remarkable as its small fine structure of 0.335 cm~^ SsJ falls 
within the range of energies displayed in Figs |2] to HI and is expected to modify the present 
conclusions. This work is currently under progress. A similar discussion obviously holds for 
the hyperfine interaction of the excited atom, which will induce even more complexity in 
the formalism. It could safely be neglected for all species except for Cesium (the hyperfine 



splitting of the 6P1/2 level is 1.167688(81) GHz (59|. 

The present formalism can be generalized to photoassociation of dipolar dimers and 
atoms, like KRb with K or KRb with Rb. If the atom is in the excited state nP and 
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the heteronuclear dimer in an excited rotational state j > 0, the long-range dimer-atom 
interaction is dominated by a dipole-quadrupole term varying as C^/R^. The long-range 
interaction between two identical dipolar molecules can also be treated in the same way, 
as the leading term will be the usual van der Waals Cq/R^ term if both molecules are in 
their lowest rotational level j = 0, or dipole-dipole C-^/R^ term if one of them is rotationally 
excited. One could thus investigate the photoassociation of two identical heteronuclear KRb 
ground state molecules in their lowest vibrational level, by a laser field with a frequency 
red-detuned with respect to j = — j = 1 transition. In this respect, photoassociation of 
two dimers is very similar to the photoassociation of two identical alkali-metal atoms, except 
that the laser frequency is much smaller for the two-dimer photoassociation. 
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